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14 . 1  INTRODUCnC^i 


Earlier  in  supersonic  aerodynamics,  it  was  stated  that  a  cotipressible, 
nonviscous,  thermally  perfect  fluid  was  adequate  for  analysis  of  problems 
(outside  of  boundary  layers)  up  to  about  Mach  =  5.0.  In  this  section  we'll 
deal  with  problems  above  M  =  5.0,  especially  as  they  relate  to  re-entry 
vehicles  and  hypersonic  cruise  vehicles. 

As  Mach  changes  from  4.99  to  5.01  nothing  dramatic  changes.  Mach  5  is 
just  a  convenient  rule  of  thunti.  Hypersonic  flew  is  best  defined  as  the 
regime  where  certain  physical  flow  phencanena  become  inportant.  The  actual 
Mach  number  Vv^iere  these  physical  pihenomena  beccane  significant  will  vary  with 
vehicle  shape,  Reynolds  number,  etc.  These  phenomena  are: 

a.  Thin,  curved  shock  layers  (i.e.  the  region  between  the  shock  and 
body) . 

b.  Strong  viscous  effects  throughout  the  shock  layer. 

c.  Lew  aensity  effects  resulting  in  a  non-zero  velocity  at  the  surface 
(velocity  slip). 

d.  High  energy  causing  the  fluid  to  behave  in  a  non- thermally  perfect 
manner. 

When  these  phenomena  become  important,  they  significantly  ccatplicate  tiie 
flow  analysis  problan.  In  addition  to  the  normal  problem  of  determining  the 
pressure  pattern  around  a  vehicle  in  order  to  get  lift,  drag  and  stability 
derivatives,  a  new  problem  becomes  iitportant  at  hypersonic  speeds;  predicting 
the  heat  input  into  the  vehicle  caused  by  the  high  energy  in  the  flow.  This 
aspect  is  not  a  ma^or  concern  at  low  Mach  numbers  but  can  be  the  dominant 
concern  at  hypersonic  speeds  ( shooting  stars  make  this  obvious ) . 

Therefore,  this  course  is  broken  into  two  broad  categories:  the 
aerodynamics  of  liigh  Mach  number  flows  (pressure  pattern,  lift,  drag)  and 
aerothermodynamics  (tanperature  predictions,  heat  transfer) .  Further,  the 
aerodynamics  of  high  Mach  flews  will  be  subdivided  into  two  sections  for 
clarity.  First,  viscosity  will  be  ignored  to  look  at  purely  Mach  number 
effects  and  tlien  the  viscous  terms  will  be  included  to  shew  their  effects  at 
lew  Reynolds  numbers. 
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14.2  MACH  NUMBER  EITECTS 


As  we  saw  in  supersonic  aerodynamics,  the  Mach  number  has  a  strong 
influence  on  shock  waves  and  the  fluid  parameter  changes  across  the  shock.  A 
hypersonic  vehicle  designed  to  achieve  orbital  velocities  will  transit  a 
region  including  Mach  numbers  up  to  at  least  M  =  25  (orbital  velocity)  and  the 
very  exi.stence  of  these  large  Mach  numbers  can  have  a  dominant  effect  on  the 
flow.  In  this  section  the  purely  fluid  dynamic  effects  of  Mach  nunber  are 
discussed  without  the  added  cortplicaticxis  of  viscosity  or  other  real  gas 
effects . 


14.2.1  Basic  Hypersonic  Shock  Relatiais 

If  we  ignore  viscosity  and  assume  a  thermally  perfect  gas  we  have  the  same 
relations  for  fluid  parameter  changes  across  a  shock  that  we  developed  in 


Chapter  6,  Supersonic  Aerodynamics.  However  for  large  Mach  nunbers,  certain 
sinplifying  assunapticns  became  possible. 

For  exaiiple,  equation  6.49  fran  supersonic  aerodynamics  gives  tlid 
relaticxiship  between  density  behind  an  oblique  shock  and  the  upstream  density: 


~  2+(Y-l)M,^sin^e 

P]  1 


(Y  t  1) 


M  1  ^  sin^fl 
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If  M  >>  1  then  M2  dcminates  so  the  2  in  the  dencminator  can  be  ignored  and 
then  the  sin^e  terms  cancel  each  other  leaving: 
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In  similar  fashion,  the  other  basic  relations  across  an  oblique  hypersonic 


shock  wave  reduce  to 
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Using  the  relationship  between  pz  and  pi  can  also  give  us  a  very  simple 
result  for  the  coeffi.cient  of  pressure,  C^.  Since  is  defined  as 

P  =  P-g— _BA  =  _2_  h  .  I 

P  qi  yMi*  Pi 


including  the  hypersonic  approximation  for  p^/Pi  gives: 
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or  with  M,  »  1, 
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In  addition,  the  relationship  between  the  wave  angle  0  and  the  wedge  angle 
6  is  sinplified.  Fran  Oiapter  6  we  had  (equation  6.55), 


tan  6 


2  cot  e 


Ml*  sin  *  6_;:;_1 _ 
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This  relationship  was  plotted  in  Figure  6.14.  For  small  wedge  angles  the 
change  in  wave  angle  gets  smaller  and  smaller  as  Mach  increases.  Figure  14.1 
is  an  illustration.  Using  small  angle  approximations: 


gives 
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if  Mj  >>  1  then  further 
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J/  *=  15 


Figure  14.1  Shock  Wave  Angle  at  High  Mach 


Qualitatively  speaking,  then,  the  soluticxi  of  hypersonic  shock  problems, 
ignoring  viscous  and  other  effects,  is  much  sinpler  than  for  supersonic  flew 
problems .  As  we '  11  see  in  the  next  sections ,  some  calculations  become 
extremely  simple  and  fortunately  give  reasonable  results  for  a  limited 
category  of  real  problems. 
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14.2.2  Mach  Number  Independence 

As  we  just  saw  in  the  preceding  section,  some  relationships  become 
independent  of  Mach  at  hypersonic  speeds.  Exaitples  are  the  density  ratio 
across  an  oblique  shock  (equation  14.1)  and  the  coefficient  of  pressure  for  a 
small  wedge  angle  (equation  14.4). 

Otlier  relationships  continue  to  depend  on  Mach  nuirber.  From  equation  14.2 
ana  14.3  it  can  be  seen  that  the  increase  in  pressure  and  tonperature  across 
an  oblique  shock  both  tend  to  »  as  M^  «>  .  If,  however,  we  non-dimensionalize 
pressure  with  dynamic  freestream  pressure. 


Pi  Vi2 


then  equation  14.3  can  be  rewritten  as 


or 
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v^hich  is  also  independent  of  Mach. 


In  general,  non-dimensional  variables  become  independent  of  Mach  at 
hypersonic  speeds.  Examples  of  this  can  be  seen  in  figures  14.2,  14.3,  and 
14.4.  In  each  of  these  exanples  the  non-dimensional  parametmeter  Cj^,  L/D^^^, 
or  C  )  all  vary  dramatically  transonically  and  supersonically,  but  caice 
hypersonic  speed  is  attained,  the  parameters  no  longer  change  with  increasing 
Mach. 

The  obvious  utility  of  this  is  ettphasized  in  the  last  example,  figure 
14.4.  Prior  to  the  first  launch  of  the  space  shuttle  no  wind  tunnel  data  was 
available  above  M  =  16  and  very  little  above  M  =  8.  While  all  parameters  do 
not  conform  to  the  predictions  as  well  as  this  example,  most  were  similar. 
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MACH 


Figure  14.2  Experimental  Drag  Measurements 

(o  Charters  and  Thomas  (1945),  x  Hodges 
(1957),  +  Stevens  (1950)) 


Figure  14.3  Wind-Tunnel  Results  for  (L/D) uay 
an  Integral  Re-entry  Body  (Krogemann,  1973) 
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Figure  14.4  Space  Shuttle  Data  for  Yawing 

Moment  due  to  Sideslip 


14.2.3  Hypersonic  Similarity 

Similarity  in  aerodynamics  refers  to  identifying  parameters  or 
combinations  of  parameters  that  result  in  identical  fluid  flew  changes  at 
different  conditions.  Ihe  advantage  is  obvious:  you  can  do  experimental 
tests  at  one  set  of  conditions  cind  then  predict  the  results  at  another  set  of 
conditions .  The  criteria  for  hypersonic  simileirity  can  be  derived 

mathematically  (see  references  2  or  3),  but  the  following  discussion  should 
suffice  to  justify  the  use  of  similarity. 

Consider  a  slender  body  in  a  hypersonic  flow  at  a  lew  angle  of  attaeik. 
These  limitations  allcw  use  of  small  perturbation  assunptions. 


SHOCK  WAVE 


6  -  WEDCE  ANGLE 


Figure  14.5  Slender  Body  in  Hypersonic  Flow 


If  we  define  the  changes  in  velocity  in  the  x  and  y  directions  as  u'  and 
v’ ,  then  at  some  point  downstream  of  the  shock, 

U  =  Voo  +  u* 

V  =  v' 

Over  a  slender  body  the  change  in  the  x  conponent  of  the  velocity  is  very 
small  ccxrpared  to  the  change  in  the  y  direction  (0  to  v'  is  a  large  percentage 
increase) .  Thus  it  can  be  seen  as  essentially  a  small  deflection  of  the  fluid 
either  up  or  dcwn.  With  this  picture. 


sin  .  =  sin  6 


But  since  u'  is  assumed  small  in  coiparison  to  V„,  then 


sin 


sin  6 
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sin  6 


or,  if  v'  «  V„ 

or 


M„  sin6_ 


Thus,  a  measure  of  the  disturbance  referenced  to  the  freestream  speed  of 
sound  is  M  sinL.  If  two  different  flow  problems  have  the  same  disturbance 
then  they  are  similar  and  have  like  solutions.  We  can  define  this  disturbance 
indicator  as  the  hypersonic  similcurity  parameter,  K. 


K  =  Moo  sin6 


(14.7) 


In  some  texts  K  is  defined  as  MooT  vhere  x  is  the  thickness  ratio,  d/jj,. 
The  thickness  ratio  for  slender  bodies  is  of  the  same  order  of  magnitude  as 
sin6  and  is  equally  good  for  stipulating  hypersonic  similarity. 
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Figure  14.6  and  14.7  are  exanples  of  data  at  different  conditions  that 
have  the  same  hypersonic  similarity  parameter.  These  examples  show  the 
usefulness  of  hypersonic  similarity  within  the  small  angle  limitations. 


Figure  14.6  Analytically  Derived  Drag  Polar  Showing 

Hypersonic  Similarity  (Ref  2:pg  50) 


PERCENT  I£NGTH 


Figure  14.7  Pressure  Ratios  for  Different  Body 

Shapes  with  Hypersonic  Similarity 
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14.2.4  Newtonian  Tneory 

In  1687,  Isaac  Newton  modeled  fluid  flew  as  a  stream  of  particles  that  do 
not  interact  with  each  other.  This  concept  is  essentially  a  totally  inviscid 
flcv/  description,  vhich  gave  very  poor  results  when  Newton  tried  to  predict 
17th  century  ship  hull  drag.  His  model,  however,  gives  surprisingly  good 
results  when  applied  to  sorte  hypersonic  problems.  Isaac  was  just  slightly 
^ead  of  his  time. 

The  basic  assunption  that  Newton  used  was  that  as  each  streamline  of 
particles  approached  a  body  the  streamline  would  be  deflected  to  parallel  the 
surface.  Essentially  this  means  that  there  will  be  a  coiplete  loss  of  normal 
nomentum  and  no  change  in  the  tangential  cerrponent  of  the  irotentum.  Figure 
14.8  shows  tills  for  a  flat  plate  at  some  angle  of  attack.  Here  V2^=  V2 

and  V,  =0. 


Figure  14,8  Newtonian  Flow  Over  a  Flat  Plate 


From  Newton's  second  law,  the  time  rate  of  change  of  norentum  is  equal  to 
the  force 


The  only  change  in  V  is  the  total  loss  of  the  nomal  oenponent  of  the 
jnomentum.  Thus  the  force  is  normal  to  the  surface  with  magnitude  : 
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F  =  m 


A  t 


For  the  flat  plate  in  figure  14.8,  this  can  be  rewritten  as 


m 


The 


F  =  — --  V,  sina 

At  X 

mass  flow  rate  incident  on  a  flat  plate  or  area  A  is 


Ttierefore 


m  =  P,  V  A  sina 


F  =  (p  V  A  sina)  V  sina 
1  1  1 


Realizing  that  i/A  is  3ust  "  P^'  we  can  rewrite  this  as 


F/A  =  p  -  P  =p  V  2  sin^a 
2  1  1 


Finally,  since  the  coefficient  of  pressure  C^,  is  just  the  static  pressure 
change  divided  by' the  dynamic  pressure,  we  get  to  Newton’s  result: 


C  =  2  sin^a 


(14.8) 


From  this  sinple  result  we  can  quickly  derive  relationships  for 
and  L/D: 


=  C  cosa  =  2  sin^a  cosa 


=  Cp  sin  =  2  sin^a 


and 


L/D  =  cycp  =  jj- 


y/hile  all  of  these  answers  gave  poor  results  for  17th  century  sailing 
ships,  figures  14.9  and  14.10  show  good  results  VN^hen  applied  to  some 
hypersonic  problems. 
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Figure  14.9.  Pressure  Distribution  on  a  Circular 

Cylinder 
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Figure  14.10  Lift-to-Drag  Ratio  for  the 

Space  Shuttle 
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Notice  that  the  non-dimensional  L/D  ratio,  Cj^,  and  Cj^  do  not  change  with 
(i/jach,  in  concert  with  the  previous  discussioi  of  Mach  number  independence. 

Newtonian  theory's  true  physical  significance  can  be  seen  ty  examining  the 
exact  oblique  shock  relationship  for  C  in  the  following  form  (see  Ref  4,  page 

Sr 

10); 


4 

y~+  I 


sin  ^6  - 


(14.9) 


This  exact  (albeit  inviscid)  equation  reduces  to  Newton's  equation  (14.8) 
v^en  Y  1  and  M  ^  -^-oo. 


14.3  VTSCDUS  hJjPLCTS 

Viscous  effects  are  ignored  whenever  possible.  Including  viscous  terms 
greatly  ccxrplicates  the  mathematical  description  of  aerodynamic  phenomena.  In 
sujiersonic  aerodynamics  we  ignored  viscosity  altogether,  relegating  it  tc  «. 
thin  boundary  larger  which  nad  a  relatively  minor  contribution  to  changes  i  i 
lift  and  drag  at  supersonic  Mach  numbers.  The  major  idea  tliat  you  should  get 
frcrn  this  section  is  that  viscous  effects  cannot  be  so  casually  dealt  with  at 
hypersonic  speeds. 

,14.3.1  Tne  Boundary  Layer 

The  boundary  layer  is  a  matliematical  sinplification.  Outside  of  some 
distance  to  the  body  in  a  fluid  flow  we  choose  to  ignore  viscosity  because  it 
sinplifies  the  flow  problem  without  causing  serious  errors  in  our  results; 
viscous  effects  are  negligible.  Inside  the  boundary  layer,  however,  viscous 
effects  must  be  included  or  our  results  will  be  significantly  in  error. 
Figure  14.11  in  an  illustration- 
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>BOmn)ABT  LAYER 


THIN  FLAT  PLATE 


Figure  14.11  Boundary  Layer 


Ihe  min  assunption  that  we  mde  in  supersonic  aerodynamics  was  that  the 
boundary  layer  thickness,  ^ ,  was  smll  conpared  to  the  scale  of  the  body. 
Hence  ignoring  i.t  caused  no  serious  errors.  However,  tlie  thickness  of  the 
boundary  layer  varies  not  only  with  distance  as  shewn  in  figure  14.11  but  also 
with  Mach  and  Reynolds  numbers  as  ; 

--  “  (14.10) 

e 

Thus  the  combinations  of  lower  Reynolds  numbers  (high  altitude  flight) 
and,  especially,  high  Mach  (thickness  grows  as  the  square  of  M»)  can  produce 
hypersonic  boundary  layers  orders  of  mgnitude  thicker  than  boundary  layers  at 
lower  speeds. 

There  are,  broadly  speaking,  two  ways  of  dealing  with  this  increased 
irtportance  of  viscosity:  use  boundary  layer  concepts  or  analytically  treat 
the  entire  flow  as  viscous.  These  are  the  subjects  of  the  next  two  brief 
sections. 
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14.3.2  ViscxDus  Interaction 


Ihe  thick  boundary  layers  at  hypersonio  speeds  lead  to  a  major  interaction 
between  the  growth  of  the  boundary  layer  and  the  outer  in  viscid  flow.  A  major 
consequence  of  tliis  interaction  is  tlie  effect  on  the  pressure  distribution,  as 
shavn  in  figure  14.12. 


X  ^ 


Figure  14,12  Pressure  Change  Due  to  Viscous  Inter¬ 
action  Over  a  Flat  Plate  at  Zero  Angle 
of  Attack 


This  increase  in  pressure  over  the  no  interaction  case  is  due  to  a  shock 
wave  created  by  the  boundary  layer  turning  and  corpressing  the  streamlines. 
If  the  flat  plate  were  at  sone  ciTigle  of  attack,  the  effect  of  the  boundary 
layer  growth  would  be  to  increase  tlie  ccrrpression  angle  and  hence  the  pressure 
increase. 

Initially,  the  grovhh  of  the  boundary  layer,  36  /3x  is  Icurge  and  thus  the 
effect  on  the  pressure  change  is  greatest.  Tnus  the  flcw'  cari  be  divided  into 
a  "strong  interaction  region"  and  a  "weak  interaction  region"  as  shewn  in 
figure  14.13.  Tire  effect  is  also  evident,  in  figure  14.12. 


Shock  Wave 


Out&r  edge 
of  b.L 


S* 


'igure  14.13  Strong  and  Weak  Viscous  Interactions 


A  detailed  analysis  of  the  viscous  interaction  phenanena  proves  that  the 
paraireter  x  is  the  governing  correlation  parameter  (ref.  3),  where 


M«  ^  w^w 

~  ~7~, — 

e 

00  oo 


(14.11) 
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From  reference  3,  for  an  insulated  flat  plate  at  zero  angle  of  attack: 

Strong  Interaction:  =  0.514X-  +  0.759  (14.12) 

Weak  Interaction:  =  1  +  0.31X  +  0.005  (14.13) 

Pco 

Experimental  verification  of  tliese  results  is  shown  in  Figure  14.14.  The 
departure  from  tlie  strong  interaction  prediction  at  low  va.lues  of  (AX)  is 
due  to  shock  separation  at  the  leading  edge  of  the  non-zero  thickness  "flat 
plate"  in  the  wind  tunnel. 


Figure  14.14  Pressure  Distribution  Induced  on  a 

Flat  Plate  at  Zero  Incidence 

Tlie  overall  effect  of  viscous  interacticxi  on  a  hypersonic  flight  vehicle 
is  to  increase  the  drag  and  hence  reduce  the  lift— to-drag  ratio. 

14.3.3  Conputational  Fluid  Dynamics 

'The  viscous  interactions  techniques  of  the  previous  section  were  developed 
in  the  1950 's  before  the  advent  of  modem  ccrputational  capability.  Ihey  are 
capable  of  providing  relatively  accurate  predictions  for  a  limited  set  of 


problems  based  on  an  approximate  theory.  The  approximation  lies  in  separate 
calculations  of  flow  characteristics  in  a  viscous  boundary  layer  and  an  outer 
inviscid  flew,  and  tben  a  coupling  of  these  separate  calculaticxis  to  take  into 
account  the  interaction  )jetween  the  two. 

The  aerodynamicist  to  calculate  viscous  effects  exactly,  sinply  by 
treating  the  entire  flow  field  between  the  body  eind  shock  as  fully  viscous, 
tto  arbitrary  division  between  a  viscous  boundary  layer  and  an  inviscid  flow 
needs  to  be  made.  These  methods  are  called  computational  fluid  dynamics 
(CFD). 

All  CFD  techniques  are  based  on  the  Navier-Stokes  equations.  The  full 
Navier-Stokes  equations  are: 


Continuity  equation 

X  momentum: 

y  momentum: 

z  momentum: 
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dr  ax  By  dz 


where 


djurj  ^  dj  UT„)  ^  djurj  ^  div7„)  divT„) 
dx  dy  dz  dx  ^  dy 

^  divT„)  djwrj  3(h>t^)  ^  djwrj 

dz  dx  By  Bz 

^  fBu  Bw\ 

r„  =  A(V.V)  +  2MT^ 
dx 

=  A(T  .  V)  + 
t,-A(V.V)  +  2m~ 

uZ 
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Wg  us  go.  extiTGiTiGly  sin^lifiGci  VGirsions  tliGSG  GQUG.tiiOii.o  in  supGirsonic 
aerod^Tiaraics .  With  CFD  fewer  and  fewer  assLinptioria  need  to  bC'  .■:ade. 

A  trulv  conplete  solution  of  the  Navi er —Stokes  equations,  however,  is 
still  beyond  our  capability.  The  irast  ccnplete  solutions  to  dale  were  done  at 

Air  Force  Flight  Dynamics  Laboratory  by  J.S.  Shane  and  S.-J.  Scneie  (ref. 
6).  They  analyzed  the  entire  flew  field  around  on  X-24C-10D  wind  tunnel  model 
at  Mach  5.95,  Re  =  1.64X10^/m,  and  a  =6°.  The  only  simplifications  to  the 
Navier-Stokes  equations  were  that  tliey  assijmed  steady-state  condizicns  so  that 
derivatives  with  respect  to  time  went  to  zero  and  they  used  a.varaged  Reynol<is 
numbers.  Even  with  these  sinplifications,  after  years  or  work  with  a  large 
ccKputer  (a  XMP-22)  they  were  only  able  to  get  solutions  for  one  steady 

state  Mach  nuntoer  at  one  angle  of  attack.  But  at  those  conditions  they  got  a 
vast  amount  of  data  which  correlated  very  well  with  the  wind  tunnel  data  tiiey 
had. 

If  you  don't  have  the  luxury  of  a  CRAY  conputer  and  years  of  time  to 
obtain  data  for  a  single  flight  condition,  CFD  is  still  useful  if  the 
Navier— Stokes  equations  are  simplified  a  little  further.  One  ccniTon 
sinplification  is  called  “parabolized  Navier-Stokes"  (PNS)  equations.  The 
full  set  of  equations  is  sinplified  by  assuming  steady— state  conditions  and 
then  dropping  the  viscous  terms  that  involve  derivatives  in  the  streanwise 
direction.  The  result  is  a  set  of  parabolic  (hence  "parabolized")  partial 
differential  equations.  These  are  much  easier  to  solve  with  a  conputer  b^'’ 
starting  witli  an  initial  data  plane  and  then  using  a  downstream  marching 
procedure.  Reference  7  gives  more  insight  into  the  use  and  utility  of  PNS 
solutions . 

Other  sinplifications  are  comonly  used  with  CFD  and  it  is  important  to 
know  the  assunptions  made  and  limitations  of  a  particular  CFD  solution.  The 
sinplest  codes  are  capable  of  very  quick  results,  the  major  difficulties  being 
input  of  the  vehicle  surface  matrix  and  extraction  of  data.  V/ith  CFD  you  get 
every  parameter  imaginable  (P,  T,  p  ,  etc.)  at  eveiry  point  in  the  flov.’. 
Interpreting  a  roemful  of  conputer  printouts  is  difficult  and  flew 
visualization  with  cenputer  generated  graphics  accessing  this  vast  data  base 
has  received  a  lot  of  attention. 
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14.4  HIGH  TI-H^b'WvlURE  EFFECl^ 

Ir.  sup-:  rt-cnic  aercx3yiTarnj.cs  \’c  deveic^x  c’  c  : 
t.aipe.r.-i.tuve  i5ccc;i  on  static  te::p€.‘rc.t:u’*e  and  a!;';irr-r.9 

cas .  V.  v’as : 


If  this  holds  true  at  hypersonic  speeds  uic;:  the  ctaca-;:.;  n  \-.fe'.perature  cn 
tlie  spe.ee  shutLle  at  M  =  25,  100,000  ft  is  51'>125  "'F*  E.y  we.y  cf  cxjnparison , 
the  surface  tanperature  of  the  sun  ic  less  than  11,000  ^F.  T3ierefore,  W3 
hopefully  conclude  that  equation  14.14  does  not  hold  ar.  iv/ix-rsonje  E]:>eeds.  2. 
say  hc^fully,  because  no  cnjrrent  materials  arc  capable  cf  vithstanding  such 
tenperatures .  Table  1  is  a.  list  of  materials  user,  in  aircraft  and  spacecraft 
construction  and  tlie  itaxiittum  tatperature  they  can  withstemd  and  still  be 
usable . 


MATERIAL 

TEMP  (°F) 

ALUMINUM 

320 

TITANIUM 

900 

INOO-  18 

1300 

RENE  41 

1600 

COLUMBIAN 

2500 

SHUTTLE  TILE 

2700 

TABLE  1.  MAXIMUM  REUSABLE  TEMPERATURE  F0f2  SEUiCr}"'')  M*\TERIALS 

The  actual  peak  tenperatures  recorded  on  the  space  rhuttie  were  under  2300 
°F,  thus  equation  14.14  and  the  assumption  cf  a  thervolly  perfect  g.as 
obviously  does  not  hold  true  at  M  “  25. 
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14.4.1  Non- Per feet  Gas 

1^10  2f0£f;.on  fCt'uetion  14.14  Goesn  t  iiolo  is  liist  et  liiglier  Macli#  snei 

hence  higher  tonrperature,  the  perfect  gas  aasu!r.ptions  break  down.  Tne  two 
main  effects  ere: 

a.  vibral'icn,  disrociation,  and  ionization  all  a’osorb  energy 
reducing  T^,,  end 

b.  The  ratio  of  specific  heats,  y  ,  decreases  with  increasing 
tearperature,  aprsrcaching  1  in  the  limi.t,  also  reducing  T^. 

For  the  first  effeert,  as  the  gas  tenperature  rises  the  force  of  impact 
between  individual  molecules  of  air  becomes  greater  and  greater.  At  sane 
point  during  tlie  increasing  taiperature,  the  force  of  individual  inpacts  is 
gr©at  enough  to  excite  internal  vi.bratcon  of  modes  of  the  molecules,  absorbing 
sane  of  the  force  (energy)  of  irrpact.  Because  of  the  normal  distriicution  Oj. 
velocities  in  a  cas,  sarv-  molecules  will  be  vibrationally  excited  at  a  lo/^er 
toiperatura,  or  velocity,  than  otliers.  Figure  14.15  shows  on  an  altitude 
velocity  map  v^here  10%  and  90%  of  air  will  be  vibrationally  excited.  As  you 
can  see,  even  in  an  SR-71  flying  at  Mach  3  and  80,000  ft  only  about  10%  of  the 
air  is  even  vibrationally  excited. 

As  we  increase  speed,  the  rising  tenperature  begins  to  tear  apart  the 
moleailar  bonds,  firsr  oxygen  and  then  nitrogen,  further  absorbing  energy. 
Much  later  indiv'idual  atoms  colliding  begin  to  knock  off  electrcxis.  Air  is 
less  than  10%  ionized  during  the  normal  space  shuttle  reentry  although  sane 
ionization  does  occur  accounting  for  the  communications  black  out  they 
experience  during  reentry'. 
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14.4.2  TemperatuT'e  ,^'5^' 

Given  the  prev.-; cs.;r-.  d;'r'-  ■  :  .v  ■  (  vh:; ;  v’c  r.c  -'.’d  very  accurately 

pin  down  the  tenTjerctur-?  c..  t.  h;  '  r  by  taking  into  acccunv. 

the  ncai-perfecrner.r  c;'  Lb;  d.'.'-c.  i-;;,?,  h~-?v;-,  a  conplicating 

i'actors.  Four  cf  are  l-  'n  y  '.'it.;  >';ynou  nur;.. 

14.4.2.1  Radiation 

bVen  with  tlie  ;voiiifyinc  eriec;.;'  yraxi cjr.ly  cef.cri'coi,  hypersonic  vehicles 
still  become  very  'no'c .  Fc:.  inrrancG,  the  Apollo  capsule  reentering  the 
atmosphere  after  retu::nir.g  frex  lun<ii-  orbit  ==  22)  hnd  peak  tenperatures  as 
high  as  11,60C  *K,  nearly  twice  the  surface  temperature  of  the  sun.  With 
taiperatures  in  tlie  fluid  this  high,  radiation  cannot  be  ignored  as  it  becemes 
the  doninant  mechanicra  for  heat  transfer.  The  equation  for  radiative  heat 
transfer  is: 

q  •  K,;T  (14.15) 

This  doesn't  seem  too  cad  until  we  realize  that  the  envLssivity  coefficient 
e ,  is  generally  a  function  cf  tenperature  also.  For  exanple,  there  are  five 
different  types  of  thermal  protection  material  on  the  space  shuttle,  all  with 
different  emissivity  coefficients.  Three  of  the  five  vary  significantly  with 
terperature  as  shown  in  Figure  14.17  (ref.  9).  The  real  mpact  of  this  is  on 
your  ability  to  predict  the  temperature.  A  small  error  in  T  can  produce  a 
large  error  in  e,  which  in  turn  v/ill  effect  T. 


14.24 


Figure  14.17  Shuttle  Tile  Emissivity  Coefficients 

14.4.2.2  Catalytic  Wall  Effects 

Another  factor  conplicating  the  prediction  of  tenperatures  is  that  near 
the  surface,  recoiibination  of  dissociated  iiolecules  or  ionized  atoms  can  occur 
at  a  higher  rate  than  in  the  equilibrium  flew.  Since  a  large  amount  of  energy 
went  into  the  dissociation  or  ionization,  it  follows  that  a  large  amount  of 
heat  may  be  released  by  the  fluid  right  next  to  the  surface  of  our  body.  The 
extent  to  v»hich  a  particular  surface  acts  as  a  catalyst  for  this  process  is 
difficult  to  predict  as  can  be  seen  in  Figure  14.18.  {ref.  10) 
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The  surface  tenperature  predicclon  assumed  a  larger  catalytic  effect  than 
was  experienced  causing  significantly  Icxver  teiiperatures  as  can  be  seen  on 
Figure  14.18  between  about  Mach  22  and  Mach  8. 

14.4.2.3  Turbulent  Boundary  Layer 

Once  the  boundary  layer  transitions  from  laminar  flcn^'  to  turbulent  the 
effect  of  conductive  heating  from  the  gas  to  the  surface  will  be  significantly 
enhanced.  Predicting  the  point  of  transition  is  very  i:rportant  to  predicting 
correct  surface  temperatures.  Again,  figure  14.18  shcv/s  that  the  transition 
is  difficult  to  predict.  The  actual  tenTperature  on  the  lower  surface  made  a 
step  increase  about  Mach  10  indicating  that  this  in  where  transition  occurred 
instead  of  about  Mach  15  vdiere  it  was  predicted. 
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14.4.2.4  Sliock  Int  eractions 


The  interacticais  of  shock  waves  off  of  various  jpai.t  '.  rf  Ik.c  vehi.cler.  cejs 
have  very  significant  effects  on  localized  heating  valuer.. 

Two  drainatic  exanples  of  this  are  the  X-  15  crd  the  ttjace  rhir  ■  le.  ITie 
X-15  was  modified  to  the  'C'  configiuration  late  in  the  pv-xr.  cjn  ella-.'ing  h.’ghej; 
Mach  numbers.  One  purpose  was  to  test  a  ramjet.  Or.  -ihc:  first  flight  of  e. 
duimy  ramjet  it  melted  off  because  of  localized  heating  cf  the  support  si.ruu 
due  to  shock  inpingement. 

The  second  example  is  the  loss  of  several  Orbital  Mcneu.;uring  System  (OMS) 
pod  tiles  on  the  early  flights  of  the  space  shuttle.  Excessive  heating 
occurred  due  to  shock  wave  interactions,  eventually  requiring  a  change  to  a 
thicker  tile  in  those  areas. 

14.4.3  Temperature  Control 

Although  accurate  prediction  of  surface  temperatures  of  hyperscnic 
vehicles  is  very  difficult,  the  designer  can  control  the  tp.nperatures  to  sore 
extent.  TVo  ways  are  ty  proper  choice  of  the  body  shape  and  the  choice  of 
reentry  profile. 

14.4.3.1  Body  Shape 

Even  though  real  air  is  not  a  perfect  gas  (fortunately),  the  vibration, 
dissociatiai,  and  ionization  that  m^e  it  non-perfect  do  not  occur 
instantaneously.  Thus  if  you  quickly  heat  up  air  (say  in  a  shock  wave)  tliei 
the  initial  temperature  will  be  very  close  to  that  predicted  assuming  a 
perfect  fluid.  Later  the  fluid  will  cool,  as  first  the  collisions  excite 
vibration  modes  and  later  as  dissociation  begins.  Figure  14.19  shews  the 
decrease. 
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In  a  moving  fluid,  the  time  axis  in  Figure  14.19  can  be  transposed  into  a 
distance  axis.  That  is,  as  the  gas  flews  past  the  shock  wave,  sane  distance 
is  required  to  allow  the  non-perfectness  of  tlie  fluid  to  absorb  sons  of  tlie 
temperature  increase  of  the  instantaneous  shock  wave.  Because  of  this,  the 
tvo  bodies  in  Figure  14.20  will  have  very  different  surface  temperatures  at 
tlie  leading  edge. 


e.  Sfaazp— nosed  b.  Blunt— xxosed 

Att&ched  Shock  Detached  Shock 

Figure  14.20  Sharp-  vs.  Blunt-  nosed  Slender  Bodies 
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The  sharp-nosed  vehiole  in  Figure  14.2Ga  will  have  surface  tenperatures  at 
the  tip  apprcaching  tVie  perfect  gas  predicti.ons.  If  the  t-iach  is  high  enough, 
say  orbital  reentry  Si?eed,  the  tip  wouid  soon  melt  and  begin  to  look  like 
Figure  14.20b. 

The  advantage  of  a  blunt-nosed  slender  body  is  that  tliS  detached  shockwave 
allows  the  air  to  internally  absorb  sane  of  the  heat  before  the  air  reaches 
the  surface  of  the  vehicle. 


14.4.3.2  Reentry  Profile 

For  orbital  reentry  vehicles,  several  choices  of  flight  path  are 
available.  Figure  14.21  s}ic.\'s  the  flight  corridor  bounded  hy  lift  and  heat 
limits  on  an  altitude-velocity  map. 

At  very  high  angles  of  attack,  range  (or  crossrange)  is  sacrificed  to 
provide  the  minirnutn  instantaneous  taiperature .  This  also  allows  you  to  expect 
a  large  tenperature  on  only  one  area  (e.g.  the  underside  of  the  shuttle  during 
a  4CP  angle  of  attcck  reentrys).  On  the  ether  haaid,  the  sliuttle's  internal 
structure  may  reach  a  higlisr  t.enperature  due  to  longer  exposure. 

Ihe  peak  temperature  occurs  just  after  landing. 


Figure  14.21  Flight  Corridor 
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